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ABSTRACT 

We present first results from a new, multiyear, time domain survey of young stars in the 
North America Nebula complex using the Palomar Transient Factory. Our survey is providing 
an unprecedented view of aperiodic variability in young stars on timescales of days to years. The 
analyzed sample covers -Rptf ~ 13.5-18 and spans a range of mid-infrared color, with larger- 
amplitude optical variables (exceeding 0.4 mag root-mean-squared) more likely to have mid- 
infrared evidence for circumstellar material. This paper characterizes infrared excess stars with 
distinct bursts above or fades below a baseline of lower-level variability, identifying 41 examples. 
The light curves exhibit a remarkable diversity of amplitudes, timescales, and morphologies, with 
a continuum of behaviors that can not be classified into distinct groups. Among the bursters, 
we identify three particularly promising sources that may represent theoretically predicted short- 
timescale accretion instabilities. Finally, we find that fading behavior is approximately twice as 
common as bursting behavior on timescales of days to years, although the bursting and fading 
duty cycle for individual objects often varies from year to year. 

Subject headings: Stars: late-type - Stars: pre-main sequence - Stars: solar-type - Stars: vari- 
ables: T Tauri, Herbig Ae/Be 



1. Introduction 

Variability in prc-main sequence stars can further our insight into physical processes associated with 
the formation and early evolution of both stars and planets. The observed flux variations are diagnos- 
tic of dynamic or radiative transfer effects that can occur on time scales ranging from hours to decades, 
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or possibly longer. Driv ing; phenomena include physical processes associated with enyelope- to-disk infall 
( Vorobvov fe Basil 2010l ), disk-to-star accretion |Wood et al. 1996 ; Mahdavi fe Kenvon 1998 ). differential 
rotat ion of a three-dimensional disk , rotation of the star, magnetic field i nteraction between th e star and the 
disk ( Romanova et al. 2004bl 20131) . accretion-driven wind and outflow ( Bans fe Konigl 2012 ). and planet- 
disk interaction. Different amplitudes and time scales can be associated with each of the theoretically 
postulated physical phenomena. In addition, the observed behavior of any individual system can be modi- 
fied by orientation with respect to the line of sight. The wide range of plausible aperiodic behavior originates 
for the most part in the circumstellar environment. Variability of circumstellar origin is superposed on an 
underlying periodic modulation that is expected due to rotation of surface inhomogeneities, analogous to 
enhanced sunspots, across the projected stellar disk, as well as any short time scale chromospheric flaring or 
other aperiodic activity. 

Historically, while many empirical studies of pre- main sequence star variability have involved searches for 
periodic behavior, those addressing aperiodic variability usually focused on explanations involving stochasti- 
cally time variable disk-to-star accretion, circumstellar extinction, or both. Examples of the former includ e 
the extreme (>2-6 mag) "outburst events" as exemplified by EX Lup and FU Ori objects (|Herbielll977|) . 
These types of sources are interpreted as undergoing episodes of rapid mass accumulation due to an insta- 
bility in the inner disk. In the context of stellar mass assembly history, the duration and frequency of such 
outbursts is important to establish empirically since these events are thought, based on theory, to play a 
determining role in setting the final mass of the star. Accretion outbursts may also deter mine a star's ap- 



pearance to us on the so-c alled "birthline" in the canonical HR diagram of stellar evolution (jHartmann et al 



1997 



Baraffe et alJl20Qgh . from which stellar masses and ages are usually derived without considering the 



effects of accretion history. For similar reasons, it is also important to understand the variati on at less ex- 



treme levels in the disk-to-star mass accretion rates. The "irregular variables" identified by, e.g. jHerbst et al 



([20021 ) at low amplitudes (<0.1-0.3 mag) and short time scales are thought to indicate non-steady accretion. 



Examples of the latter, extinction-related, variability include UX Ori stars, which undergo distinct and 
somewhat long-lived extinction events, as well as the broader category of stars identified by color-color 
and color-magnitude trends consi stent with shorter-timescale, random variation along reddening vectors hy 



Morales- Calderon et al 



e.g-.ICarpenter et al.l (I2001L 120021) ■ More recently, so-called "dipper" events (e.g.. lCodv fe Hillenbrandll201 



20111 ) are attributed to repeated sub-day or several-day circumstellar extinction en- 



hancem ents. Such repeating but ape riod ic flux dips or eclipse-l ike events have been qualitatively explained 
by, e.g.. lFlahertv fc Muzerollel (|2010h and lFlahertv et al.l |2012bl) using rotating non-axisymmetric disk mod- 
els or bv lTurner et al.l ([20101 ) with a vertical disk turbulence model. Periodic versions of the dipper class are 
known as AA Tau stars (e.g., Bouvier et al. 19991) . 



While variability is a common property of young pre-main-sequence stars and has been viewed as a key 
observational characteristic of newly formed stars since their discovery (e.g.. |Jovlll945[ ). the full breadth of 



variable phenomena has not been explored in quantitative detail. The expected periodic variabil i ty, indicative 



■P 1 

of ste l lar rotation, has been w ell studied fe.g. jGrininl l2000l:lHerbst et al.ll2007l: llrwin et al.ll2008l; iRebull et al 



2006c ICieza fc Baliberl 120071 ) . The harder to interpret aperiodic variability, our focus, is cataloged but 
relatively unexplored in the literature. Many fundamental properties of aperiodic variability in young stars 
are still poorly quantified. As aperiodic variables constitute more than half to two-thirds of variable stars in 
star-forming regions, characterizing and understanding them is essential to completing our understanding of 
young star/disk physics. In addi tion to the photometric stu dy reported here , spec troscopic monitoring may 
be required, such as reported bv lChoudhurv et al.l (|201ll ) or ICostigan et al.l ([20121 ). 



To further progress, we have carried out a multiyear optical wavelength monitoring program aimed at 
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determining the distributions of variability time scales and amplitudes among members of a young stellar 
populatio n. Specifically, we have su rveyed several square degrees of the North America and Pelican Nebulae 
region ([Reipurth &: Schneider 2008 ) using the wide field of view and reliable time coverage of the Palomar 
Transient Factor^ (PTF: lLaw et all boifj iRau et alJliooj) in operation at the Palomar 48" Samuel Oschin 
Telescope. A notable niche of these PTF data is the long duration of the time series at roughly nightly 
cadence. Future contributions will assess the overall variability statistics, categories of variability, and 
characteristic time scales and amplitudes of optical variability in the North America and Pelican Nebula 
region. Here, we present our PTF survey strategy and our methodology for identifying variable stars. We 
then investigate two specific types of variability phenomena exhibited among an infrared-excess selected 
sample of objects. 

The present paper focuses on observable optical variabilit y among the ~ 210 known and suspected 
members of the North America Nebula complex cataloged by iRebull et al.1 (|2011l ) based on mid-infrared 
selection techniques. Of these, 84% are within our monitored field. Among the wide range of behaviors ex- 
hibited by variable stars, we consider the evidence for and typical properties of bursting or fading behavior, 
possibly mixed with other forms of variability. In the case of bursting stars, while accretion- related insta- 
bilities having timescales of a few tens of day s have been predicted by a number of theoretical studies (e.g., 
Alv fc Kuiipenlll990 ; iRomanova et al. |2004aL no evidence for accretion bursts prod uced by such instab ilities 
has bee n published (Bouvier et al.112007 ). although accretion bursts on both shorter ( Murphy et aljEoil ) and 
longer (|HerbisJl2008l ) timescales has been observed. We assess the frequency of these intermediate timescale 
instabilities. For fading stars, while the existence of extinction-related variability is well-established, results 
vary among authors as to the frequency of young stars exhibiting such behavior, as well as the typical 
timescales. We also address in this study the ratio of bursting to fading light curves for a typical T Tauri 
star population. 

The remainder of this paper is organized as follows. In Section [21 we present our photometric data and 
our detection thresholds for variability. Section [3] discusses how we defined the burster and fader populations 
and their key properties. In Sections |4] and El we discuss the bursters and faders in more detail, with an 
emphasis on how the largest sample yet identified of such objects can constrain their underlying physics. In 
Section[6l we describe several noteworthy objects in more detail. Section [7] summarizes our results, describes 
limitations of our analysis, and suggests pathways for future work. 



2. Photometric Data 

2.1. Instrument and Survey 

The PTF survey camera is a mosaic of 11 chips, covering a total area of 7.26 square degrees with 1" 
pixels. This is a wide enough field to observe almost the entire North America Nebula complex in a single 
exposure, as shown in Figure [TJ PTF typically observes in the Mould- R or g' bands. Because our targets 
are intrinsically red and further reddened by extinction, we conducted our entire survey in R band, where 
a typical 60-second exposure reaches stars as faint as 20th magnitude. Throughout the survey, we took at 
least two exposures per night, separated by one hour. 

Our survey cadence was complex as a result of changing operational factors. Our observations started 
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in August 2009, continuing with observations every third night until October, when Palomar was shut down 
due to ash from fires. When we started our 2010 season in April, the cadence was lowered to every fifth 
night. From August to October 2010, we were able to observe every night, while the remainder of the season 
was hampered by poor weather. For our 2011 and 2012 seasons, from March 2011 to January 2012 and from 
March 2012 to December 2012, we were able to observe every night, but only during bright time, as the 
PTF project had started observing exclusively with the g' filter during dark time. In addition, during July 
and August 2011, we obtained hourly exposures all night, in both bright and dark time. We illustrate our 
observing pattern in Figure [2j 

Our survey represents one of the most uninterrupted, multiyear optical variability surveys of a star- 
forming region, featuring 884 epochs across 352 nights between 2009 and 2012. Our largest data gaps are the 
1-3 month gaps in the winter, when the region is not visible from Palomar, and the two-week gaps during 
dark time in most months when the R filter is not available. Aside from these regular gaps, we have two 
years of uninterrupted nightly coverage, except for occasional weather gaps, and another year and a half of 
lower-cadence data for probing long-term variability. 



2.2. Reduction Pipeline 



All the PTF data for our field was processed by the PTF Photometric Pipeline. Images were debi- 
ased, flatfielded, and astrometrically calibrated, with source catalogs generated by SExtractor (Laher et al., 
in prep). An absolute photometric calibratio n good to a syst ematic limit of ~ 2% was generated using 
SDSS sources observed throughout the night (jOfek et al.l 120121 ). Relative photometric calibration further 
refi ned the pho t omet ry, pa rticularly on n o nphot ometric nights (Levitan et al., in prep; for algorithm details 



see lOfek et al.l (|2011l ) and iLevitan et al.1 (|201ll )). The PTF Photometric Pipeline photometry is typically 
repeatable to 0.5-1% for bright (15th mag) nonvariable sources in sparse fields on photometric nights. Pho- 
tometry for typical sources in our field is less reliable, of the order of 2-3%, because nebula emission and 
source crowding introduce additional errors. 

The pipeline flagged photometric points as bad detections if the sources were automatically identified 
as part of airplane, satellite, or cosmic-ray tracks; if they fell on a high dark current, unusually noisy, or 
poorly illuminated area of the chip; if they fell on a chip edge; if they contained dead pixels; if they were 
affected by bleeding from bright stars; if they contained saturated pixels; or if they had neighbors biasing 
their photometry. We removed any sources from our sample that were flagged in more than half the epochs, 
as discussed in the next section, and we removed all flagged epochs from our light curves before plotting or 
analyzing them. 



2.3. Identifying the Variables 

To determine which sources were variable during the observation period, we grouped all PTF detections 
with 14 < R < 20 mag into half-magnitude bins on a chip-by-chip basis. The width of the bins (0.5 mag) 
was chosen so that the brightest and least populated bin (14-14.5 mag) had roughly 100 sources on most 
chips. We then computed the median RMS of all the stars in each bin, and fit the medians by an equation 
of the form: 

RMS = J a 2 + (bx 10-0.4p(m O9 -14))2 (!) 
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This equation is partly motivated as the sum of a systematic term and a flux-dependent term; the exponent 
of the flux-dependence p was allowed to vary because the natural choice, V — \ (i- e -i noise that scales as the 
square root of the flux, as expected from photon noise), was too shallow. In practice we found p ~ | for 
most chips. We list the fit parameters in Tabled] 

The curve found by fitting Equation [T] describes the locus of nonvariable stars on a given chip. We 
defined the boundary between variable and nonvariable stars to be 1.75 times the median RMS. This cutoff 
was determined empirically, rather than analytically, to avoid making assumptions about the noise properties 
of the data. We set the cutoff by visually inspecting light curves with both R ~ 14 and R ~ 16; at RMS 
values lower than 1.75 times the threshold the light curves were indistinguishable from noise, while at higher 
values the light curves were clearly structured on short timescales. 

We show in Figure|3]plots of RMS vs. magnitude for the six chips that covered the star-forming complex, 
along with the median fit and the variability detection boundary for each chip. For 14th magnitude stars, 
we are sensitive to variability with an RMS amplitude of a few percent, while below 16th magnitude, we can 
probe only 10% flux variations. In Table [TJ we list the number of PTF sources and the number and fraction 
identified as photometrically variable using the methods outlined above. Nearly 3000 variables projected on 
the dark cloud and the associated nebulae are identified. Their RMS amplitudes range from 0.03 to 1.1 mag. 



3. Bursting and Fading Among Infrared Excess Sources 
3.1. Sample Selection 



Because the North America Nebula complex is located in the plane of the Galaxy, a significant number of 
our high quality light curves are those of foreground or background field stars. In our first reconnaissance of 
the variability properties of the region, we ther efore concent r ated o n variable sta rs among a lis t of candidate 
North America Nebula members identified by iRebull et al.1 (|2011l ) . Specifically, iRebull et al.1 used infrared 
colors, primarily Spitzcr IRAC 3.6 /im — MIPS 24 /im, to identify stars surrounded by circumstellar dust. 
Additional considerations included location in various color-magnitude diagrams that help distinguish young 
stars from contaminating dusty sources such as extragalactic AGN and galactic late-type giants. Each 
source was assigned a spectral energy distribution class based on the slope of a linear fit to all available 
photometry between 2 and 24 /im. Class I sources have rising slopes and are interpreted as objects with 
not only circumstellar disks, but likely more spherically distributed envelopes as well. Flat-spectrum sources 
have roughly constant \F\ over the 2-24 fim. range and have a similar interpretation. Class II sources are 
consistent with traditional disk SEDs. Class III sources have the steepest slopes; most have no excess in the 



IRAC bands but were selected based on an excess at 24 fim. Only 6 of the Class III sources in IRebull et al 



([201 11 ) were not selected using either IRAC or MIPS excess criteria. IRebull et al.1 note that, since their 
primary selection is based on infrared data, they are incomplete with respect to Class III sources. 



Of the 2082 candidates from IRebull et al.l ([201 lh . 601 had a counterpart in the PTF source catalog. As 
we show in Table [2] the recovery rate by PTF depended strongly on the type of IR excess. Only 5% of the 
relatively red Class I sources in the PTF field had detections, while fully 93% of the relatively blue Class III 
sources were detected by PTF. The strong correlation between (infrared) source color and recovery rate, in 
the sense that redder sources are recovered less often, suggests that most of the sources we did not recover 
in PTF were missed because they were below our optical detection limits. However, we also know, from 
image inspection, that the PTF pipeline had difficulty identifying and extracting sources from crowded or 
nebulous regions. If we assume that all the Class III sources must be bright enough to detect in the optical 
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if they are visible in the Spitzer bands even with a small infrared excess, then source extraction problems 
should dominate the 7% missing Class III sources. Presumably, roughly 7% of the rest of the sample also 
fell in regions where the PTF pipeline could not reliably identify sources. We note that, while the overall 
incompleteness does not affect our main science goals, the bias away from Class III sources in the parent 
sample and the bias away from Class I sources from cross-matching to PTF do limit our ability to examine 
how variability properties change with the degree of infrared excess. 

From our sample of 601 infrared excess selected candidate members with PTF counterparts, we restricted 
our attention to the 253 sources brighter than a median i?pTF = 18. The detailed breakdown by SED type is 
given in Table [2j We found from experience that the photometric quality for sources fainter than -Rptf ~ 18 
was such that, while we could determine whether a source was variable, we could not consistently assess 
the structure of the variability. Considering only sources whose light curves had bad photometry flags (see 
Section l2~2l for a list) in fewer than half the epochs further reduced the sample to 186 stars, which are shown 
in Figure HI The figure shows no trend with -Rptf except for more sources at fainter magnitudes, suggesting 
our magnitude limits avoid any substantial systematics. High-amplitude sources (RMS > 0.3-0.4 mag) tend 
to be associated with strong infrared excess, while low amplitudes are found in both strong- and weak-excess 
sources. 

From this sample of 186, we studied in more detail the 117 that showed significant variability in PTF, 
as defined in Section 12.31 Both cuts are presented in more detail in Table [2j These infrared excess selected 
variables include most of the high amplitude variables in the field, as shown in Figure [5j most of the low 
amplitude variables in the field lack an infrared excess and are not part of our sample. The 117 infrared- 
excess variable sources, along with the other variables in the field, exhibit a wide range of time scales and 
amplitudes in their light curves. We sought to categorize the lightcurves and hereafter we focus on those 
that can be identified as bursting or fading. 

When selecting sources for inclusion on the list of bursters or fadcrs, we defined a burst in a light curve 
as a period of elevated fluxes above a (local) floor of relatively constant brightness. We did not place any 
explicit restriction on the length of the candidate burst. However, we tended to require elevated fluxes in 
multiple consecutive epochs to be certain that a brighter measurement was not a measurement error, and 
we required that the period of elevated fluxes be short enough that we could recognize the remainder of the 
light curve as a well-defined "quiescent" state. We defined fades analogously: a period of lowered fluxes, 
with the caveats that we believed the lower fluxes represented real variability and that the lower fluxes were 
distinct from the normal variability of the star. Both definitions were necessarily subjective, and we review 
possible selection effects in Section [731 

We visually inspected all 117 light curves for bursting or fading activity. For comparison, we also 
inspected 100 randomly chosen variable PTF sources that did not have an infrared excess, mixing them with 
the sample of 117 so that we did not know whether any particular light curve was from the target sample 
or the control group. We designated a star as a burster or a fader if it had at least one bursting or fading 
event during the monitoring period. 

3.2. Burster and Fader Statistics 

We identified 14 stars with candidate bursts and 29 stars with candidate fades, with two stars showing 
both bursting and fading behavior. The sources are listed in Table [3l with their photometric behavior 
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summarized in Table HJ Light curves of all 41 stars are available online from the PTF websitcQ The sources 
are also highlighted in Figures [TJ |31 and @] 

For comparison, in the control group of 100 sources with no infrared excess, we saw only two stars that 
appeared to have one burst each, and no faders other than eclipsing binaries. The burst detected in one of 
the stars turned out to be a transient scattered light artifact we had failed to spot at the time of the original 
analysis. The other may also have been identified as a burster because of a systematic error in the data or 
in our visual inspection, or it may represent real astrophysical variability in the field. In the former case we 
expect ~ 2 of the bursters in our target sample to be mislabeled, while in the latter we expect ~ 1 false 
positives. 

The stars listed in Tables [3] and SJ some of which are highlighted in Figures [6] and El show a wide variety 
of behaviors. We see variability from a few tenths of a magnitude to several magnitudes. The bursts or 
fades last anywhere from around a day, the shortest timescale resolvable in most of our data, to hundreds of 
days. Events may repeat as frequently as once a week, or can appear only once in the three- year monitoring 
period. Nearly all the bursters and faders are aperiodic, with the exception of two faders that are discussed 
further in Section [SJ 



3.3. Spectroscopic Characterization 



We pursued optical spectroscopy of both the variable star selected sample (this paper) and the infrared- 
excess selected sample of iRebull et al.l (|2011[ ) using the MMT, Keck Observatory, Palomar Observatory, and 
Kitt Peak National Observatory. 

We observed 164 vari able infrared-exce ss sources in the North America Nebula using the DEIMOS 
multi-object spectrog raph (jFaber et al.ll2003l) at Keck on 2012 July 18-19, using the 600 line/mm grating. 
The sample included 19 bursters or faders. PTF was monitoring the field during both nights that spectra 
were taken, allowing us to determine the photometric state represented by each spectrum for all stars except 
those varying significantly in less than a day. 

The spectra were reduced using a modified version of the DEEP2 pipeline ([Newman et al. boii l ICooper et al 



2012J), provided to us by Evan Kirby. The spectra were bias-corrected, dome-flatficlded, and lamp-calibrated, 
but were not flux-calibrated. We corrected for sky and nebula emission by subtracting the off-source spectrum 
visible within each slit. The final spectra covered approximately the 4400-9500 A range at 5 A resolution, 
although the range covered by the spectrum of any particular star could shift by ~ 500 A in either direction 
depending on the position of the star's slit on the instrument mask. Many cosmic rays were left uncorrected 
by the pipeline, so when making the figures in this paper, we cleaned the cosmic rays by hand for clarity. 

194 sources sele cted by either variabi lity or infrared excess were observed using the Hectospec multi- 
object spectrograph (jFabricant et al.l 120051) on the MMT on 2012 November 3, December 4, and December 
6, using the 270 lines/mm grating. The sample included 22 b ursters or faders. The data were pipeline 
processed at the Harvard-Smithsonian Center for Astrophysics ([Mink et al.1 120071 ). The final spectra cover 
3700-9100 A at 6 A resolution. PTF observed the region on November 3 and December 6, but to interpret 
the December 4 spectra we had to interpolate between photometry from December 3 and December 6. 



One of the authors (L. H.) had previously obtained low-resolution optical spectra of sources in the 



2 http : //www. astro . caltech.edu/ptf/ 
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North America Nebula complex with the HYDRA multi-object spectrograph (jBarden et al.lll993l ) on the 
3.5m WIYN telescope at Kitt Peak, using the 316 line/mm grating, on six nights between 1998 June 2 
and 1998 July 21. L. H . also took spectra using the (now decommissioned) Norris multi-object spectrograph 
([Hamilton et al.lll993l ) on the 5m Hale Telescope at Palomar, using the 600 line/mm grating, on 1998 August 
14-15, 1999 July 21-23, and 1999 September 2-5. We took spectra of 27 bursters or faders during these runs. 
The HYDRA and Norris spectra do not have concurrent photometry. 

The HYDRA and Norris observations were reduced for us using custom routines written in IDL. The 
routines applied corrections for bias, dome flats, cosmic rays, scattered light, and wavelength calibration. The 
spectra were not flux-calibrated. Sky and nebula emission were corrected by taking a shorter sky exposure 
offset 6-10" from the target position, and subtracting the counts in the sky exposure from the corresponding 
source spectrum, after scaling to the difference in observing times. In several configurations the sky emission 
was scaled by an additional 10-20% to account for changes in the sky transmission. The HYDRA spectra 
covered 5000-10000 A at R ~ 1500, while the Norris spectra covered 6100-8750 A at R ~ 2000. 



4. The Burster Phenomenon 

4.1. Population Properties 

Upward excursions in young star light curves traditionally have been assumed to come from one of two 
mechanisms. Long-lasting, several-magnitude events in young stars with circumstellar accretion disks (e.g., 
EX Lup, FU Ori) are interpreted as dramatic increases in the accretion rate from the disk to the star. Events 
lasting a few hours or less and rising by at most a few tenths of a magnitude, particularly in disk-free stars, 
have been assumed to be associated with magnetic flares like those seen on the young field star UV Cet or 
on the Sun. White light flares on the Sun last tens o f minutes, while those on M dwarfs last up to several 
hours (e.g., Kowalski et al. 20ld Kretzschmar 2011 ). As these timescales are set by the cooling times of 



dense chromospheric material at the base of the coronal loop, it is unlikely that magnetic flares can produce 
optical emission with much longer durations than observed. 

Of the 14 stars that show bursting behavior, only two, [OSP2002] BRC 31 8 and FHO 1, have bursts 
lasting 1-2 hours, short enough to be plausible flares. The remainder must be driven by temporary increases 
in accretion, drops in extinction, or some other phenomenon. The bursters show a wide variety of behaviors. 
None are strictly periodic, although [OSP2002] BRC 31 8 and FHO 29 do show enhanced photometric activity 
at roughly 300-day intervals. Some bursters, like FHO 26, repeat every few weeks. Others, like LkHa 185 
or FHO 4, show bursts only once a year or even more rarely. While [OSP2002] BRC 31 8 and FHO 1 have 
very short bursts, too brief to resolve outside our highest-cadence monitoring in mid-2011, FHO 17 featured 
a burst lasting over 100 days, and FHO 18 showed bursts with a range of lengths from a few days to two 
weeks. 

Despite their variety, the bursters do not fall naturally into distinct subclasses, forming instead a contin- 
uum of behaviors. We show in Figure [8] the joint distributions of burst amplitudes, burst widths, and burst 
separations for all 14 bursters. To avoid systematics associated with separating a burst or fade from the 
surrounding, sometimes complex, variability, and to avoid complications from varying sampling from event 
to event, the timescales and amplitudes in Figure [5] were estimated by eye and should be taken as illustrative 
values only. There is no pattern visible in the plot aside from a rough trend where longer bursts tend to 
be separated by longer intervals. The absence of distinct groups of bursters suggests that the diversity of 
sources can be explained by continuously varying the parameters of a single common scenario, rather than 
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by invoking different mechanisms or different configurations for short- and long-timescale bursters. 



If either enhanced accretion or reduced extinction are responsible for bursting events, then stars with 
large infrared excess, and therefore more circumstellar material, may be more likely to show bursting behavior 
than stars with small infrared excess. Using the Kendall's r statistic ( Kendall 1938f ), we found no evidence 
for a correlation between the Spitzer IRAC/MIPS [3.6] — [24] color (56% confidence) and the presence of 
bursting among our sample of 117 IR-excess sources, but we did find a marginally significant correlation 
(2.4% confidence) with Spitzer IRAC [3.6] — [8.0] color, in the sense that redder sources are more likely 
to show bursting behavior. We note that of the 14 bursters, only two, FHO 2 and FHO 24, are Class III 
sources. The rest have K — [8] > 1.8 and K — [24] > 5 (see Figure|4]for comparison to the rest of the sample). 
We note that since the K-band and Spitzer data are not coeval, the reported colors may be distorted by 
variability between the epochs of observation. However, mid-infrared variability is typically a few tenths of a 
magnitude or less (jEspaillat et al.ll2011t Morales- Calderon et al.ll2011t iFlahertv et al.ll2012al ). and so should 
not dramatically affect a star's position in diagrams such as Figure 2) 

While the weak correlation with [3.6] — [8.0] color suggests that bursters are associated with stronger 
circumstellar disks, and therefore with the possibility of enhanced accretion or reduced circumstellar extinc- 
tion, the absence of a similar correlation with [3.6] — [24] color weakens this result. As noted in Section |3~T1 
however, only a limited range of infrared color is well-represented in this sample. It is also possible that any 
correlation is being diluted by radiative transfer effects, geometry, or other factors that determine whether 
any particular star shows bursting behavior. We discuss how additional data could allow more conclusive 
tests in Section [7731 



4.2. Constraints on Short-Term Accretion Outbursts 



Magnetic or viscous instabilities acting at the boundary between the stellar magnetosphere and the 
circumstellar disk are expected to prod uce short bursts of accretion on timescales of weeks to months for 
certain regim es of disk properties (e.g., lAlv fc Kuiipersl Il990t iGoodson fc Winded Il999l: iRomanova et al 



2004al . l2005h . However, variability from such outbursts has never been observed (jBouvier et alJ 120071 ). The 
consistent cadence and long time coverage of our PTF survey have allowed the most sensitive search to date 
for such accretion events. 

Of the bursting sources in Table FHO 2, FHO 4, and FHO 24 show multiple bursts lasting tens of 
days each. The separations between bursts vary: tens of days in the case of FHO 24, 100-300 days in the 
case of FHO 2, and ~ 500 days for FHO 4. We show all three sources in Figure The timescales and 
shapes of t hese events, particularly FHO 2 and FHO 4, resemble the simulated variations in M shown in 
Figure 4 of IRomanova et alJ (|2004aI K Although they do not stand out in the context of our sample, where 
burst durations vary continuously from < 1-150 days, FHO 2, FHO 4, and FHO 24 are noteworthy as the 
first bursts reported in young stars having timescales of tens of days. To our knowledge, these light curves 
represent the first observations consistent with the predicted inner-disk instabilities. 

Models predict that short-term accretion outbursts should have amplitudes of a few tenths of a mag- 
nitude. For example, scaling to a fiducial star with 0.8 M© and 2 R@, the simulations of IRomanova et al 
(|2004al ) predict an accretion rate of 2 x 10 -8 M^yr " 1 in quiescence and 6-8 x 10 -8 M Q yr _1 in outburst. The 



fiducial star would have a luminosity of 1.14 L© (Sie ss et al ] l2000h . with quiescent and outburst accretion 
luminosities of 0.25 L© and 0.88 L©, respectively, implying a brightening of ~ 0.4 magnitudes between 
quiescence and outburst. The three candidate stars have amplitudes between 0.2 and 0.5 mag, consistent 
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values given star-to-star variations in star and disk parameters. 

The behavior of the light curves is inconsistent with white-light flares analogous to those seen on the 
Sun. White-light flares tend to have a steep rise followed by an exponential decay. None of the three 
bursts show an exponential profile, and the timescales of tens of days are much longer than the minutes 
to hours durations observed in the Sun or in low-mass stars. A superposition of many short flares is also 
unlikely: FHO 2 and FHO 24 show little variability on timescales of a single day, as might be expected from 
a stochastic sum of shorter events. In addition, a 0.4 mag burst lasting 30 days corresponds to an energy 
release, depending on the (unknown) spectrum of the transient emission, of ~ 6 x 10 39 erg for a 1.4Lq star. 
The entire stellar magnetic field, integrating a dipole field from an assumed i?* ~ 2Rq to infinite radius, 
contai ns only ~ 5 x 10 38 ( B BUTi /l kG) erg. Even a 3 kG field, near the largest values observed in T Tauri 



stars (|Bouvier et al.ll2007l ). cannot provide enough energy to power the bursts. 



The brightness enhancements for these three objects are also unlikely to be dust-clearing episodes. Only 
FHO 4 is a Class II source, while FHO 2 and FHO 24 are both Class III sources with excess only in the MIPS 
24 /xm band. It is doubtful that these two stars have enough circumstellar dust in the inner disk to allow 
significant extinction-driven variability. We note that stars with infrar ed excess only in the MIPS bands can 
still s how ongoing accretion on the order of 0.1-0.5 x 10~ 8 Af Q yr _1 ( Muzerolle et al. 2009t Espaillat et al 



2012I ). so the absence of an IRAC excess does not rule out either low- level accretion or, plausibly, brief 
periods of accretion at a higher rate. Extinction is a possible origin for the variability of FHO 4; color data 
taken over the course of one of its bursts could test this hypothesis. 



The Fader Phenomenon 



The t wo prototypical faders are AA Tau, which fades repeatedly by 1.4 mag over 30% of an 8.2- 
day cycle ( Bouvier et al. 19991 2003 ). and UX Ori, which fades by 3 mag for tens of days at irregular 



intervals jWaters fc WaelkenJl998h . Both AA Tau and UX Ori are well understood as the result of recurring 
extinction by circumstellar material, from a warped inner disk edge in the case of AA Tau or from more 
irregular structures in the case of UX Ori. 

Of the 29 sources that show some kind of fading behavior, only two, LkHa 174 and FHO 12, show the 
periodic modulation characteristic of AA Tau. Four more sources, LkHa 150, FHO 7, FHO 15, and FHO 27, 
show multiple fading events with durations of tens of days, as seen for UX Ori stars. However, their typical 
amplitudes of 1 mag or less are much smaller than the 3-4 mag fades associated with UX Ori stars. 

The remaining 23 faders do not resemble either of the previously established categories. The natural 
assumption is that these sources also have their variability dominated by circumstellar extinction, with 
different spatial scales or different geometries causing the light curves to behave differently. 

All the sources except for the two AA Tau analogs are aperiodic, and, as illustrated in Figure they 
often bear little resemblance to each other. For example, FHO 19 has narrow fades repeating every 8-10 
days, but without enough coherence to be periodic. In contrast, NSV 25414 and FHO 3 both have frequent 
but irregular events, with the interval between adjacent fades varying by more than a factor of two. At 
the other extreme, FHO 21 and FHO 22 each show only one fading event per year, while LkHa 150 and 
FHO 25 fade only once in the entire survey period. Most fading events are short, but those of LkHa 150 and 
[OSP2002] BRC 31 1 last for hundreds of days. Most stars have fading events of roughly constant depth, but 
FHO 15 and FHO 20 have significant amplitude variability. Most fading events are symmetric, but FHO 11 
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and FHO 27 show strongly lopsided events. 

Like the bursters, the faders do not separate naturally into sources with distinct timescales. We show 
in Figure |S] the joint distributions of fade amplitudes, fade widths, and fade separations for all 29 faders. 
The absence of gaps in the plot suggests that, as with the bursters, short- and long-timescale faders have a 
common origin. 

As with the bursters, we tested for a correlation with the presence of circumstellar material. Using 
the Kendall's r statistic, we found no evidence for a correlation between the Spitzer IRAC [3.6] — [8.0] or 
IRAC/MIPS [3.6] — [24] colors and the presence of fading among our sample of 117 IR-excess variables, at 
20% and 18% confidence respectively. However, we do not find a single example of a fader among Class III 
sources (i.e., significant excess only at 24 fj.m), as would be expected if circumstellar material near the star 
is needed for fading events to occur. 

To test whether the fading events could instead be the result of variable foreground extinction, we 
searched for a correlation between stars' near-infrared color, where we can avoid variability-induced system- 
atic errors through the use of coeval 2MASS photometry, and the presence or absence of fading behavior. 
Since unreddened M dwarfs have J — K < 1, stars with 1 < J — K < 3 must have significant extinction, 
while stars with J — K < 1 may have only moderate extinction. If fading events are caused by foreground 
dust, we might expect fading to be more prevalent among the reddest stars. Using the Kendall's r statistic, 
we found no evidence for a correlation between the J — K color and the presence of fading among our sample 
of 117 IR-excess variables, at 26% confidence. 

While we find that neither the degree of infrared excess nor proxies for near-infrared reddening are good 
predictors for the presence of fading behavior, the absence of faders among the Class III sources is consistent 
with fading events requiring the presence of inner disk dust and therefore the possibility of occasionally 
enhanced extinction along the line of sight. We discuss how additional data could allow more conclusive 
tests in Section [731 

6. Individual Sources of Interest 

In Sections |4] and [5] we have examined the 41 burster or fader candidates as an ensemble. However, 
many of the sources have a character of their own. While we present brief descriptions of all the sources in 
Table |4j in this section we focus on a small number of stars whose behavior seems particularly difficult to 
explain. We present the available data on each and challenge interested readers to develop models for these 
sources. 



6.1. FHO 26 

FHO 26 showed several-day-long, ~ 0.7 mag bursts in 2010 and 2011 (see upper right panel of Figure [T0|) 
but became quiescent in late 2011. In 2012, except for two brief bursts, it has shown only a 0.2 mag, 5.6-day 
periodic modulation. The 2010-2011 bursts do not phase up under the 2012 period. FHO 26 has a modest 
infrared excess, as shown in the lower left panel of Figure HUl 

We show in the lower right panel a spectrum of the source taken in July 2012, well into the quiescent 
phase and near the peak of the periodic variability. The spectrum shows an M4.5 photosphere with emission 
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from Ha (-13 A equivalent width). 



6.2. [OSP2002] BRC 31 1 

[OSP2002] BRC 31 1 grew fainter by nearly three magnitudes between April and August 2011 but showed 
relatively little variability before the fade, as shown in the upper left panel of Figure [TT] Our spectrum, 
taken during the star's faint state, shows a forest of emission lines including Ha, Ca II, [O I], [Fe II], [S II], 
[Ni II], Fe II, and many others. A spectrum of [OSP2002] BRC 31 1 from 1998 shows only Ha, Ca II, and 
Fe II at the same strength as in 2012, plus much weaker [O I] and [Fe II] lines. We see few absorption lines 
in the spectrum in either epoch. 

The 1998 and 2012 spectra are similar to high- and low-state spectra, respectively, of the long-term 



variable PTFlOnvg (|Hillenbrand et al.ll2012j ). Since BRC 31 1, like PTFlOnvg, is a Class I infrared excess 
source, it is possible that BRC 31 1 is a similar system: a high-inclination source with circumstellar material 
obscuring the inner disk, stellar photosphere, and accretion zone, b ut not obscur i ng a s patially extended jet. 



We note the light curve resembles that of VI 184 Tau presented bv lGrinin et al.l (|2008l ). 



6.3. FHO 18 

FHO 18, shown in the upper panels of Figure [T^l faded twice by 0.6 mag in quick succession in 2011. 
Immediately before each fading event, it brightened by 0.3 mag (upper right panel). This behavior was 
not repeated for other fading events during our monitoring period. Aside from these two fades and their 
precursor bursts, FHO 18 appears to be a typical Class II young star. 

Our DEIMOS spectrum of FHO 18 was taken during a 0.4 mag fade. The spectrum shows a K5 star 
with Ha emission (-23 A equivalent width) as well as weaker Ca II and He I emission. However, as this 
fade was not preceded by a burst the spectrum does not directly constrain the star's unusual behavior in 
mid-2011. 



6.4. FHO 27 

FHO 27 had only 0.5-0.6 mag variability with a roughly constant or slightly rising mean magnitude 
throughout 2009-2010, but then began to show deep (up to 2 mag) fading events from late 2011 onward. At 
the same time, the upper envelope of the light curve began to gradually dim, leveling off in mid-2012 after a 
total decrease of roughly 0.8-0.9 mag. The minimum magnitude reached during each fading event rose from 
17.9 in the first fade to 17.3 in late 2012, so that the most recent fades have only a few tenths of a magnitude 
depth. 

While the fading events repeat every 30-90 days, they are not periodic. In addition, each fade has a 
dramatically different profile from those before and after, with many of the profiles showing strong asym- 
metries (Figure H3l upper right panel), and some fades being much shallower than the rest (e.g., a 0.3 mag 
fade in late April 2012 was sandwiched between a 1.1 mag and a 0.6 mag fade). 

The spectr al energy distributio n of the source, shown in the lower left panel of Figure [13l has a strong 



infrared excess; iRebull et al.l (|2011l ) classify FHO 27 as a fiat-spectrum source. 
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We acquired one spectrum of FHO 27 in July 2012, during the star's long-term low state but between 
the deeper fading events. The K7 spectrum in the lower right panel of Figure [T31 shows strong Ha (-80 A 
equivalent width), Paschen series, and Ca II emission. Weaker lines in the spectrum include He I, [O I], and 
O I. 



6.5. FHO 28 

Like FHO 27, FHO 28 was dominated by 0.6 mag irregular variability in the first few years of our survey, 
interrupted by occasional 1 mag fades. Then, in early 2012, it began showing rapid variability with the same 
maximum brightness, but with a much higher amplitude of 2 mag. The high-amplitude variability lasted 
130 days before the star returned to its earlier behavior. Since the source was not strictly periodic during 
its high-amplitude phase, it is not clear whether the variability has been fully resolved at our daily cadence, 
in which case the fades are roughly 9 days apart, or whether we are seeing a strobing effect of a more rapid 
23-hour variation. The light curve is shown on the top two panels of Figure 1141 

The spectral energy distribution, shown in the lower left panel of Figure [Ml shows a Class II infrared 
excess. A spectrum of FHO 28 (Figure H~4l lower right panel), taken during its strongly varying phase, shows 
an M3 star with strong Ha emission (-60 A equivalent width) and weak Ca II lines. An older spectrum 
shows that Ha was much weaker (-20 A) in 1998, although since we don't know the photometric state of 
FHO 28 at the time, it is not clear whether the difference between the two spectra is related to the star's 
increased activity in 2012. 

FHO 28 is yet another example of how the photometric behavior of young stars can change abruptly 
from one year to the next. This source would not be classified as a fader if we had only data from its 
active phase, as the photometry shows no preference between high and low magnitudes (Figure 1141 upper 
right panel). It is the comparison to previous years that allows us to establish that the brighter magnitudes 
represent the unperturbed state of the star. 

7. Summary and Discussion 
7.1. Key Results 

We have presented first results from a new variability survey of young stars that probes a large dynamic 
range of timescales, from roughly a day to roughly a year. We have used this new data set to uniformly identify 
stars with episodic variability, regardless of whether the episodes had day-scale, month-scale, or year-scale 
durations and regardless of whether the episodes were periodic. From a sample of 186 candidate members of 
the North America Nebula complex, we have identified 14 that showed episodes of brightening ("bursters") 
and 29 that showed episodes of dimming ( "faders" ) , sometimes mixed with more erratic variability. Two stars 
showed both bursting and fading at separate epochs. We have presented basic photometric and spectroscopic 
properties of both bursters and faders. 

We have found that: 

1. Most high-amplitude variables have a strong infrared excess, while low-amplitude variables may or may 
not have a strong excess. While similar correlations have been noted before, i.e., that classical T Tauri 
stars tend to have higher amplitudes than weak-lined T Tauri stars, we show here that a correlation 
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between degree of infrared excess and variability amplitude also holds among stars with infrared excess. 

2. Even within the individual burster and fader classes, we see a wide variety of timescales, amplitudes, and 
burst or fade profiles. This includes events that occur only once or twice in our three-year monitoring 
period, and would be missed in a shorter survey. It is not clear whether these varied behaviors imply 
varied underlying mechanisms. We find no gap separating groups of bursters or faders with different 
amplitudes or timescales (Figure [8]), suggesting that they are all members of a single population, but 
in-depth study of representative objects will be needed to settle the issue. 

3. We identify three bursters whose photometric and spectroscopic characteristics are consistent with 
published models of accretion driven by instabilities at the boundary between the stellar magnetosphere 
and the circumstellar disk. To our knowledge, this is the first time candidate objects corresponding to 
these models have been identified. 

4. A substantial number of sources show variability over long timescales. Among other examples, FHO 14 
and FHO 28 showed enhanced fading activity in an interval 100-200 days long. [OSP2002] BRC 31 1 
changed from a 15th magnitude star in 2010 to a 18-19th magnitude star in 2012. [OSP2002] BRC 31 8 
and FHO 29 both showed bursting modulated by a timescale of roughly 300 days. Except for the sudden 
decay of [OSP2002] BRC 31 1, these are behaviors that have not been associated with bursting or fading 
activity before, for lack of sufficient sampling. 



7.2. Comparison to Previous Work 



While much previous time domain work on young stars has focused on finding and characterizing periodic 
variables, there have been some studies of more general variability. Here we discuss whether our population 
statistics are consistent with the existing literature. 

We see bursting behavior in 14 sources, or 12 ± 3% of the R < 18 mag variables with an infrared excess. 
To the best of our knowledge, no one has reported a short-term optical burst fraction for pre-main-sequence 
stars, so we have no published results to which to compare this figure. 

We see fading behavior in 29 sources, or 25 ± 4% of the variables with an infrare d excess and 1 6 ± 3% 
of all sources with a good PTF light curve and an infrared excess. For comparison, Ulencar et all |2010h 
found 28% of stars sele cted from X-ray or Ha emissi on, all variables, showed periodic fading behavior in 
unfiltcred optical light. iMorales-Calderon et al.l (|201ll ) found fades (periodic or not) in the mid-infrared in 
5% of variables and 3% of their total sampl e, selected by proper motion, X-ray or Ha emission, or infrared 
excess. Finally, ICodv fc Hillenbrand! (l2010h found I-band fading behavior in 6% of their variables and 5% 
of their total sample, selected by kinematics, Ha emission, forbidden line emission, lithium absorption, or 
infrared excess. Each of these surveys was a few weeks in duration, shorter than our survey, but had higher 
cadence by factors of 10-200. 

To test whether our results are consistent with previous work after accounting for differences in our ob- 
serving strategies, we clipped our light curves to a 30-day period of high-cadence obs ervations, up to eight per 



night, between JD 2455765.5 and 2455795.5. This allowed us to compare our data to lMorales-Calderon et al 



((20111), who observed their field for a month at roughly a 2-hour cadence. We found that 12 of our faders 
(LkHa 174, V1701 Cyg, and FHO 3, 5, 14, 16, 18, 19, 20, 21, 22, and 28) were recognizable as such during 
the 30-day period, indicating that with only a month of high-cadence data we would have reported a 10 ± 3% 
fader fraction out of the variables in our sample or 6 ± 2% of the infrared-selected sample. This is slightly 
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higher than, though consistent with, the iMorales-Ca ldcron ct al. results. Si nce our ground-based survey 



had more data gaps than the Spitzer observations of iMorales-Calderon et all however, our fader rate had 
we observed with their exact cadence may have been higher. On the other hand, it is possible that we are 
overestimating our recovery rate, since we had already identified these stars as faders using the full data set 
and were aware of their nature while examining the clipped light curves, introducing hindsight bias. 



7.3. Limitations of the Present Work 

We were careful to identify bursting and fading events using only the light curves themselves, and not 
any ancillary data such as SEDs or spectra, to avoid psychological biases in interpreting ambiguous cases. 
However, we could not eliminate all ambiguity: the qualitative nature of event identification inevitably made 
some kinds of events easier to detect than others. The easiest events to identify were either those where 
the event lasted for several days, so that the light curve resolved the event profile, or those where the event 
repeated many times over our observing baseline, so that we could be confident that a high or low point 
represented real variability rather than a statistical fluke or an isolated error in the data reduction. We tried 
to confirm, using thumbnail images, whether isolated high or low points represented brief but real changes 
in the stellar flux, but image inspection allowed us to verify only high-amplitude events. We therefore may 
be incomplete to variability on timescales of a few hours. 

We also had difficulty identifying bursts or fades lasting longer than several months, particularly if they 
were superimposed on other variability. Some stars in our sample showed erratic variability on timescales of 
months, and it is not always clear from only three years of data whether a star that spent several months in 
a high (or low) state had undergone an anomalous change in brightness, or whether we were merely seeing 
an extreme in a continuous series of brightness fluctuations. We chose to err on the side of caution, and only 
counted sources where the light curve apart from the candidate burst or fade had much lower-amplitude 
variability. However, this introduced a bias against mixed variability modes. 

There are at least two sources in the North Am erica Nebula comp lex that, while they meet our definition 
of bursters, are absent from our sample. PTFlOqpf ( Miller et al.ll2011 ) was an R — 16.5 star at the beginning 



of the survey that brightened to R ~ 12.5 in mid-2010 and has remained there since. The source was 
disqualified from this paper's sample because it failed three criteria in the photometry produced by the 
PTF Photometric Pipeline: it was flagged as blended with nearby stars at nearly all epochs, it was flagged 
as saturated in nearly a ll epochs after the outburst, and its median magnitude of 12.9 was well above our 



flux limits. PTFlOnvg (jCovev et al.ll2011l) did not rise past PTF's saturation limits; however, as noted in 
Section I3TT1 the PTF Photometric Pipeline had difficulty identifying sources around nebulosity. PTFlOnvg 
is located just off a bright nebula filament, and neither it nor any other nearby sources were extracted. 
These two omissions illustrate key sources of incompleteness in this work: crowding, nebula contamination, 
and a limited magnitude range. Fortunately, these problems do not apply to the majority of sources in the 
Spitzer-selected sample, which are well-separated, in low-background regions, and of less than one magnitude 
amplitude. 

This work is based primarily on a long-term, single-band photometric survey, which has allowed us to 
identify and characterize new types of bursters and faders. However, the data we have presented here cannot 
identify without ambiguity the physics behind each kind of bursting or fading variability, or even whether all 
bursters or all faders represent different cases of a common variability mechanism. The following additional 
data would provide more insight into the nature of bursters or faders: 
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• Time series color information would help test whether all of the faders are caused by variable extinction 
along the line of sight to the star. Color data would also help us interpret bursters by providing color 
constraints to better estimate the energy released in the burst. We plan to present a color analysis of 
our bursters and faders in future work. 

• Spectroscopic monitoring, especially at high dispersion, would allow us to compare accretion and wind 
indicators in a star's high and low states, allowing us to distinguish which events are accretion-powered, 
which represent partial or total obscuration of the photosphere, and which are driven by something 
else entirely. 

• Polarimetry would help identify which bursts or fades are associated with changes in the obscuration 
of the star, as it probes what fraction of the measu r ed flux comes from the photosphere and what 



fraction is scattered from the disk (e.g., Grinin 1992 ; Bouvier et al. 19991 ). In particular, it could be 
used to directly test the hypothesis that all fades are obscuration events - if they are, then they should 
all show stronger polarization at minimum light. 



We have shown that the class of faders is far broader than previously appreciated, and that bursters, 
while fewer in number, show a comparable diversity. We have identified new phenomenology within both 
classes. These objects can serve as prototypes for future study of particular forms of bursting or fading 
activity. 
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Fig. 1. — The North America Nebula complex, as observed by PTF in a single epoch from 2009. Only six 
of the 11 PTF chips are shown; the remainder, to the left of this field, were off the nebula and probed the 
Galactic field population. The North America Nebula proper (NGC 7000) is on the left side of the frame, 
while the Pelican Nebula (IC 5070) is on the upper right, with the L935 dark cloud between them. Th e 
blue circles mark the positions of candidate members selected using infrared excess by iRebull et al.l ([201 11 ) . 
We also highlight stars with apparent bursting activity from Table [4] with red X's, and stars with apparent 
fading activity with red squares. 
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Fig. 2. — The complex cadence used in the PTF North America Nebula survey. Date labels denote the 
beginning of each year. Left: histogram of the number of observations taken in each week of the survey. 
Right: individual observation times of all survey epochs, with each month dispersed along the vertical axis 
for clarity. The period of all-night, high-cadence monitoring appears as a set of elongated points in mid-2011. 
The cadence in 2011-2012 was close to nightly, while observations in 2009 and 2010 were more sporadic. 
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Fig. 3. — RMS scatter vs. median magnitude for all sources with flags (listed in Section |2~2|) in fewer than 
half the epochs. The fit to the median RMS as a function of magnitude is plotted as the lower red curve, 
while our variability detection threshold (1.75 times the median) is plotted as the upper red curve. As in 
Figure [TJ X's mark candidate bursting stars while squares mark candidate fading stars. 27 high-amplitude 
variables are beyond the upper edge of the Chip plot, 10 each above the upper edge of the Chip 1 and 
Chip 2 plots, and 1-5 off the upper edge of each of the others. 
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Fig. 4. — PTF mag nitude and IR color distributions for those PTF sources that have an infrared excess from 



Rebull et alJ (| 2 1 lh and whose light curves have flags (listed in Section I2.2[) in fewer than half the epochs. 



The color of the dots indicates the degree of infrared excess: blue dots are class III sources, green ones 
class II, yellow ones have a flat IR spectrum, while magenta sources are class I sources. The black sources 
are those that were not detected in the Spitzer 24/mi band, and so did not have an IR excess class listed in 



Rebull et al.1 (|201lf ). Not all sources appear on both plots, as som e had missi n g 8^tm or 24/im photometry. 



The curves in the upper left panel show synthetic photometry of ISiess et all (pOOOh isochrones for ages of 
2 Myr (red) and 100 Myr (blue), at a distance of 600 pc, indicating the expected colors of stars with no 
infrared excess at all. As in Figures [T] and [3l X's mark candidate bursting stars while squares mark candidate 
fading stars. 
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Fig. 5. — RMS distribution of the sample, and its correlation with the presence of Spitzer infrared excess. 
Left: the fraction of sources with an infrared excess out of all PTF sources with 13.5 < R < 18 and flags 
(listed in Section I2.2|) in fewer than half the epoc hs. Right: the solid line denotes the distribution of RMS 
amplitudes for the 186 infrared excess sources from lRebull et al.l ([20111 ) whose light curves have flags in fewer 
than half the epochs. The dashed line represents the subset of 117 that lie above the variability thresholds 
in Figure [H 
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Fig. 6. — Examples of the diverse behavior seen in our lightcurves for bursting stars. Variability amplitudes 
range from a few tenths of a magnitude to 2 magnitudes (top row). Detected bursts can last from less than 
two days to over a hundred (middle row), and can be separated by anywhere from 10 days to a year (bottom 
row). For scale, the horizontal bar near the top of each panel shows a 10 day interval. No points having any 
of the flags listed in Section l2~2l are plotted. 
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Fig. 7. — Examples of the diverse behavior seen in our lightcurves for fadcrs. Variability amplitudes range 
from a few tenths of a magnitude to nearly 2 magnitudes (top row). Fades can last anywhere from one day 
to over a year (middle row), and can be separated by anywhere from 9 days to over a year (bottom row). 
For scale, the horizontal bar near the top of each panel shows a 10 day interval. No points having any of the 
flags listed in Section l2~2l are plotted. 
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Burst Width (days) Fade Width (days) 

Fig. 8. — Plots of the characteristic amplitudes and timescales for bursters (left) and faders (right), illustrat- 
ing the wide variety of observed events. The plus signs indicate sources that have either well-defined fixed 
values for their amplitudes and timescales (such as the two periodic AA Tau analogs in the lower left corner 
of the fader panels), or that have only a single measurement (such as the single-event sources in the upper 
right of any panel). The ellipses represent sources that have bursts or fades of varying amplitudes, varying 
widths, or varying separations within a single light curve. The area below the dotted line on the lower two 
figures, where events would need to overlap each other, is not allowed, though some ellipses appear there 
because this analysis does not consider correlations between width and separation. 
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Fig. 9. — Burst profiles for three bursters with durations of tens of days. No points having any of the flags 
listed in Section [2.2l are plotted. In this and subsequent plots, the points are connected by line segments to 
clarify the order of closely spaced observations. FHO 2 has the simplest event profile, showing a smooth rise 
and fall over a 20-30 day interval. The bursts of FHO 4 are longer and show a more complex profile. The 
light curve for FHO 24 shows a large number of contiguous bursts rather than a few isolated events like the 
other two. 
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Fig. 10. — A star whose regular bursting activity stopped at the end of 2011. The upper left panel shows the 
full 3- year light curve, with the shaded period expanded in the upper right panel to illustrate typical bursts 
for this source and the vertical red line marking the time at which the July 2012 spectrum in the lower right 
panel was taken. The red scale bar represents a 10-day interval. No points having any of the flags listed 
in Section |2"T21 are plotted. The lower left panel shows the spectral energy distribution for this source. The 
points are taken from non-simultaneou s optical, near-infrared , and Spitzer photometry. The solid curve is a 
reddened NextGen model atmosphere ( Hauschildt et al. 19991) with temperature corresponding to the star's 
spectral type, matched to the optical through J-band fluxes. 
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Fig. 11. — Same as Figure [101 but for a star showing a sudden decline in 2011. The star was not detected 
in roughly half the epochs in 2012; the non-detections are not shown. The upper right panel highlights the 
decline, including a temporary dip that interrupted it. Since we could not determine a spectral type for this 
source, the photosphere shown in the lower left panel is for an assumed effective temperature. 
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Fig. 12. — Same as FigureUHl but for a star showing an odd combination of bursting and fading. The upper 
right panel shows the fades and their precursor bursts. 
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Fig. 13. — Same as Figure HU1 but for a star showing a series of fades superimposed on a year-long decay. 
The upper right panel highlights the asymmetric profile of one of the fades. 
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Fig. 14. — Same as Figure \TU\ but for a star showing an increased frequency of fades in 2012. The upper 
right panel illustrates that the fades in the more active phase were nearly superimposed. 
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Table 1. Fits to RMS as a Function of Magnitude by Chip 



Chip 


a 


6 


P 


Total Sources 


Variables 





0.033 


0.0053 


0.68 


6090 


491 (8.1%) 


1 


0.018 


0.0130 


0.52 


9330 


490 (5.3%) 


2 


0.023 


0.0145 


0.49 


14057 


533 (3.8%) 


G 


0.026 


0.0072 


0.62 


11036 


457 (4.1%) 


7 


0.025 


0.0081 


0.60 


7720 


337 (4.4%) 


8 


0.022 


0.0105 


0.55 


10393 


397 (3.8%) 



Note. — Best- fit values for the parameters in Equation [T] for 
each of the six target chips. The parameter a can be interpreted 
as the systematic noise floor at bright magnitudes, and p is the 
power-law dependence of RMS on flux at faint magnitudes. The 
last two columns show the total number of PTF sources on each 
chip as well as the number selected by making an RMS cut at 
1.75 times the value given by Equation [T] 



Table 2. Selection Effects in the PTF Sample 



IR Source Type 


# in RGS2011 


in PTF Field 


PTF Counterparts 


With R < 18 


Flags in < 50% Epochs 


High RMS 


MlPS-Only 


25 


25 (100%) 


( 0%) 








Class I 


273 


242 ( 89%) 


11 ( 5%) 


3 (27%) 


3 (100%) 


3 (100%) 


Flat-Spectrum 


272 


242 ( 89%) 


53 (22%) 


20 (38%) 


13 ( 65%) 


10 ( 77%) 


Class II 


604 


542 ( 90%) 


321 (59%) 


160 (50%) 


120 ( 75%) 


79 ( 66%) 


Class III 


112 


82 ( 73%) 


76 (93%) 


43 (53%) 


31 ( 72%) 


16 ( 52%) 


IRAC-Only 


796 


613 ( 77%) 


140 (23%) 


27 (19%) 


19 ( 70%) 


9 ( 47%) 


Total 


2082 


1746 ( 84%) 


601 (34%) 


253 (42%) 


186 ( 74%) 


117 ( 63%) 



Note. — iRebull et~ai] l l201ll) gave SED classes only for sources that were detected in both IRAC and MIPS. Sources that were detected 
in only one or the other are listed for comparison, but were not used to estimate the incompleteness from source confusion and flux limits. 
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Tabic 3. Statistics of the Candidate Bursters and Faders 



[RGS2011] ID 


Short 


SED 


Ft 


Rmed 


RMS 


AR 


Total 


Unflagged 




Name 


Class 


(mag) 


(mag) 


(mag) 


(mag) 


Detections 


Detections 


205032.32+442617.4 


FHO 1 


II 


17 A 


17.5 


0.13 


1.12 


873 


870 


205036.93+442140.8 


[OSP2002] BRC 31 1 I 


16.8 


16.6 


1.24 


4.94 


750 


505 


205040.29+443049.0 


LkHo 139 


II 


14.6 


14.6 


0.13 


0.92 


883 


778 


205042.78+442155.8 


[OSP2002] BRC 31 8 


16.9 


16.9 


0.18 


1.97 


877 


872 


205100.90+443149.8 


V1701 Cyg 


II 


15.5 


15.4 


0.36 


1.89 


769 


629 


205114.80+424819.8 


FHO 2 


III 


13.8 


13.8 


0.075 


0.91 


857 


750 


205115.14+441817.4 


LkHa 150 


II 


16.4 


16.3 


0.29 


2.22 


879 


833 


205119.43+441930.5 


FHO 3 


11 


17.1 


17.0 


0.48 


2.74 


873 


865 


205120.99+442619.6 


LkHct 153 


II 


15.1 


15.1 


0.13 


0.80 


882 


866 


205123.59+441542.5 


FHO 4 


II 


17.6 


17.6 


0.13 


1.13 


875 


875 


205124.70+441818.5 


FHO 5 


II 


18.0 


17.9 


0.24 


2.06 


872 


839 


205139.26+442428.0 


FHO 6 


II 


14.8 


14.7 


0.26 


1.49 


884 


872 


205139.93+443314. 1 


FHO 7 


11 


16.7 


16.5 


0.53 


2.66 


879 


877 


205145.99+442835.1 


FHO 8 


II 


17.8 


17.7 


0.27 


2.71 


874 


870 


205155.70+443352.6 


FHO 9 


II 


15.9 


15.8 


0.26 


1.76 


881 


865 


205158.63+441456.7 


FHO 10 


Flat 


16.7 


16.7 


0.12 


0.90 


879 


825 


205203.65+442838.1 


FHO 11 


II 


18.0 


17.9 


0.13 


0.98 


870 


836 


205228.33+442114.7 


FHO 12 


II 


16.5 


16.4 


0.26 


1.71 


878 


863 


205230.89+442011.3 


LkHa 174 


II 


16.7 


16.7 


0.28 


1.51 


878 


878 


205252.48+441424.9 


FHO 13 


II 


18.1 


18. l a 


0.42 


2.63 


874 


795 


205253.43+441936.3 


FHO 14 


II 


18.0 


17.9 


0.14 


1.11 


873 


853 


205254.30+435216.3 


FHO 15 


II 


17.1 


17.0 


0.34 


2.01 


877 


785 


205314.00+441257.8 


FHO 16 


II 


17.1 


17.0 


0.18 


1.21 


877 


877 


205315.62+434422.8 


FHO 17 


II 


17.3 


17.5 


0.49 


2.21 


848 


720 


205340.13+441045.6 


FHO 18 


11 


17.0 


17.0 


0.21 


1.87 


875 


875 


205410.15+443103.0 


FHO 19 




18.0 


17.9 


0.24 


1.88 


869 


867 


205413.74+442432.4 


FHO 20 


II 


16.3 


16.2 


0.21 


2.36 


876 


876 


205424.41+444817.3 


FHO 21 


II 


16.7 


16.6 


0.24 


1.59 


876 


876 


205445.66+444341.8 


FHO 22 




17.4 


17.3 


0.31 


2.99 


874 


872 


205446.61+441205.7 


FHO 23 


II 


17.4 


17.3 


0.16 


0.87 


763 


665 


205451.27+430622.6 


FHO 24 


III 


15.9 


15.8 


0.13 


0.70 


860 


860 


205503.01+441051.9 


FHO 25 


Flat 


16.1 


16.1 


0.12 


1.31 


876 


872 


205534.30+432637.1 


[CP2005] 17 


II 


17.2 


17.2 


0.10 


0.88 


850 


850 


205659.32+434752.9 


FHO 26 




18.0 


18. a 


0.23 


1.72 


861 


827 


205759.84+435326.5 


LkHo 185 


11 


14.6 


14.6 


0.074 


0.70 


884 


884 


205801.36+434520.5 


FHO 27 


Flat 


16.6 


16.3 


0.61 


2.35 


878 


866 


205806.10+435301.4 


V1716 Cyg 


II 


16.5 


16.5 


0.16 


1.10 


879 


879 


205825.55+435328.6 


FHO 28 


II 


17.8 


17.7 


0.40 


2.40 


874 


871 


205839.73+440132.8 


FHO 29 


Flat 


16.7 


16.8 


0.59 


3.17 


879 


877 


205905.98+442655.9 


NSV 25414 


II 


14.7 


14.6 


0.45 


2.21 


884 


823 


205906.69+441823.7 


FHO 30 


II 


17.2 


17.2 


0.14 


1.23 


872 


869 



Note. — R denotes the mean PTF magnitude, R me d the median PTF magnitude, and AR the peak-to-pcak amplitude. 

a While this star is fainter than R me d = 18 in the latest data release, the target selection was done using an earlier release, 
at which time the source had R me d < 18. 
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Table 4. Phenomenology of Candidate Bursters and Faders 



[RGS2011] ID 


Short 




Event 


Lightcurve 


Spectrum 




Name 


(mag) 




Notes 


Notes 


205032.32+442617.4 


FHO 1 


17.5 


Burster 


Several bursts, each lasting 












only 1-2 hours. 




205036.93+442140.8 


[OSP2002] 


16.3 


Fader 


Faded in mid-2011, still at min- 


Spectrum dominated by emis- 




BRC 31 1 






imum. See section [6.21 


sion lines in both 1998 and 












2012. Both epochs show Ha, 












Ca II, Paschen series, and Fe II 












[,(ii (i i i r i i xi 

emission; 2012 also has [O 1], 












[Fe II], [S II], and [Ni II]. 


205040.29+443049.0 


LkHa 139 


14.6 


Burster 


One burst lasting 3 days in 












2011. 




205042.78+442155.8 


[OSP2002] 


16.9 


Burster 


300-day modulation, with 






BRC 31 8 






daily 1-2 hour bursts near 












maximum of the modulation. 




205100.90+443149.8 


V1701 Cyg 


15.3 


Fader 


Fades lasting several days, 












roughly once a month. 




205114.80+424819.8 


FHO 2 


13.8 


Burster 


Bursts lasting ~ 50 days, every 












100-300 days. See Section [4751 




205115.14+441817.4 


LkHa 150 


16.3 


Fader 


Faded by 1 mag in early 2012 












for 3-4 months. Long rise with 












±0.4 mag variations during re- 












covery. 




205119.43+441930.5 


FHO 3 


16.9 


Fader 


2-day fades at intervals from 4 












to 7 days. 




205 120. 99+4426 19.0 


LkHa 153 


15.0 


Burster 


One burst lasting 2-15 days, 












and several lasting less than 1 












day each 




205123.59+441542.5 


FHO 4 


17.6 


Burster 


Two bursts lasting ~ 60 days, 


M2 star with Ha, He I, [N II], 










separated by 350 days. More 


Ca II emission. Ha and Ca II 










complex profile than FHO 2. 


half as strong in 2012 as in 










See Section [4721 


1998. 


205124.70+441818.5 


FHO 5 


17.9 


Fader 


Many short 1 mag fades last- 












ing ~ 1 day, mixed with some 












longer (~ 3 day) but shallower 












(~ 0.6 mag) fades. 




205139.26+442428.0 


FHO 6 


14.8 


Fader 


Many short fades lasting ~ 












4 days, separated by 20-50 












days, superimposed on lower- 












amplitude erratic variability. 




205139.93+443314.1 


FHO 7 


16.3 


Fader 


Three fading events, 50, 150, 


K5 star with Ha and Ca II 



and > 180 days long, each with emission, 
complex bursts in their cores; 
one additional event, lasting < 
70 days, started at the end of 
the 2010 season. Decline seen 
over 2012, though not as pro- 
nounced as in FHO 27. 
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Table 4 — Continued 



[RGS2011] ID 


Short 


H-med 


Event 


Lightcurve 


Spectrum 




Name 


(mag) 




Notes 


Notes 



205145.99+442835.1 FHO 8 
205155.70+443352.6 FHO 9 



205158.63+441456.7 FHO 10 



205203.65+442838.1 FHO 11 



205228.33+442114.7 FHO 12 



205230.89+442011.3 LkHa 174 



205252.48+441424.9 FHO 13 



205253.43+441936.3 FHO 14 



205254.30+435216.3 FHO 15 



17.7 Fader 

15.9 Fader 

16.7 Fader 

Burster 

17.9 Fader 

16.5 Fader 

16.7 Fader 

18.0 Fader 

18.0 Fader 

17.1 Fader 



205314.00+441257.8 FHO 16 



17.1 Fader 



One 0.5 mag fading event last- 
ing 100 days. 

One 0.6 mag fading event 
lasting ~ 6 days in July 
2011. 4-day, 0.6 mag, fading 
events separated by 10-20 days 
throughout rest of light curve. 
Two fades ~ 0.3 mag deep last- 
ing 4-5 days, and one lasting < 
10 days. Events separated by 
several months. Mixed with er- 
ratic variability of ~ 0.2 mag. 
Two bursts ~ 0.4 mag high 
lasting < 3 days each, sep- 
arated by 7 days. Mixed 
with erratic variability of ~ 
0.2 mag. 

Slow decay over ~ 100 days 
followed by a rapid rise in ~ 
30 days. Weaker, shorter fade 
2 years before had a fast decay 
followed by a slow rise. 
1.5-day fading events repeating 
every 5.8 days. 

Fading events lasting 3 days, 
repeating every 7.7 days. 
Roughly 1/3 of the cycles do 
not have a fade. 
Fades lasting several days, ev- 
ery 10-20 days. Most fades 
have depths of ~ 1 mag; 
roughly every ~ 200 days a 
fade is deeper, ~ 1.4 mag. 
0.4-0.7 mag fading events last- 
ing 6-12 days every 20-30 days. 
One 0.2 mag, 150-day fad- 
ing event with several of the 
shorter fades within it. 
Three low states lasting 100, 
30, and 70 days, in order, sep- 
arated by one year; all three 
show high variability at mini- 
mum. First two fades 1.3 mag 
deep, third only 0.8 mag. 
Combination of 0.6 mag fades, 
lasting 2-4 days, and 0.3 mag 
fades, lasting 60-80 days. 



K7 star with strong Ha emis- 
sion, as well as He I and [O I] 
emission. 

K5 star with Ho, Ca II, and 
He I emission. 



K8 with Ha, He I, and O I 
emission. 
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Table 4 — Continued 



[RGS2011] ID 


Short 




Event 


Lightcurve 


Spectrum 




Name 


(mag) 




Notes 


Notes 



205315.62+434422.8 FHO 17 



205340.13+441045.6 FHO 18 



205446.61+441205.7 FHO 23 



205451.27+430622.6 FHO 24 



205503.01+441051.9 FHO 25 



205801.36+434520.5 FHO 27 



17.6 



17.0 



205410.15+443103.0 FHO 19 18.0 

205413.74+442432.4 FHO 20 16.2 

205424.41+444817.3 FHO 21 16.6 

205445.66+444341.8 FHO 22 17.3 



17.3 



15.9 



16.0 



205534.30+432637.1 [CP2005] 17 17.1 
205659.32+434752.9 FHO 26 17.9 



205759.84+435326.5 LkHo 185 14.1 



16.1 



Burster 
Fader 

Burster 

Fader 

Fader 
Fader 

Fader 



Fader 



Burster 



Fader 

Fader 
Burster 



Burster 



Fader 



Several 0.4 mag bursts lasting 

1- 3 days, followed by a qui- 
escent period, followed by a 
1.5 mag burst lasting 150 days. 
Two 0.4 mag fades lasting 5 
and 3 days, 11 days apart. 
Both fades immediately pre- 
ceded by 0.3 mag bursts. See 
section 16.31 

Two 0.8 mag bursts lasting 10 
and 7 days, 240 days apart. 
Several 0.3 mag bursts sepa- 
rated by tens of days. 
Several fades lasting 3 days 
each, repeating every 8-10 
days. Fade depth varies be- 
tween 0.5 and 0.9 mag. 

2- 5 day fading events; longer 
events tend to be deeper. 
Three fades, lasting ~ 10 days 
(first part not observed), 5 
days, and 11 days, separated 
by 250 and 330 days. 
Complex fades lasting 6-20 
days, separated by 230 and 300 
days. Hints of a double pro- 
file for each event. One addi- 
tional 3-day fade 50 days after 
the third main fade. 

Several fades lasting 2-6 days, 
separated by a few weeks. 
Fades range from 0.6 mag to 
0.3 mag, the level of the under- 
lying erratic variability. 
0.2 mag burst lasting ~ 15 
days in 2010, followed by a se- 
ries of 0.5 mag bursts in 2012 
lasting 15-40 days each. See 
Section l4~2l 

One ~ 5-10-day fade in late 
2010 

One 65-day fade in 2010. 
Several bursts in 2010-2011, 
lasting 4-5 days each and sepa- 
rated by 10-30 days. No activ- 
ity in 2012. See section IB~T1 
First half of a 0.3 mag burst 
before a data gap in mid-2011. 
Rise time 2 days. 
Multiple fading events lasting 
15-40 days and separated by 
intervals ranging from 30-60 
days. Events superimposed on 
a steep decline over the course 
of 2012, more extreme than in 
FHO 7. Fading events get shal- 
lower over the course of the de- 
cline. See section [6. II 



K5 star with Ha, He I, Ca II, 
and [N II] emission. 



M4.5 star with Ho emission. 



K7 star with strong Ho, Ca II, 
Paschen series, and He I emis- 
sion, and weaker lines of [O I] 
and O I. 
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Table 4 — Continued 



[RGS2011] ID 


Short 


D 


Event 


Lightcurve 


Spectrum 




Name 


(mag) 




Notes 


Notes 


205806.10+435301.4 


V1716 Cyg 


16.5 


Burster 


Two bursts, the first lasting 5- 












20 days and the second 3 days, 












separated by 35 days. Com- 












plex profiles. 




205825.55+435328.6 


FHO 28 


17.7 


Fader 


130-day interval of repeated 


M3 star with Ha emission in 










8-day fading events in 2012; 


both 1998 and 2012, though 










only 5, well-separated events 


the line is stronger in 2012. 










each in 2010 and 2011. 2011 


The 2012 spectrum also has 










fades were typically only 2 days 


weak emission of Ca II, [N II], 










long, while 2010 events were 


He I. 










too sparsely sampled to con- 












strain their length. See section 












1531 




205839.73+440132.8 


FHO 29 


16.8 


Burster 


High states in early 2010, early 












2011, late 2011, and entire first 












half of 2012. 2010-2011 bursts 












repeat roughly every 270-300 












days, but 2012 behavior does 












not fit the period. 




205905.98+442655.9 


NSV 25414 


14.6 


Fader 


1 mag fading events lasting 10- 












15 days, with ±0.5 mag vari- 












ability at minimum. Fades re- 












peat every ~ 30 days. 




205906.69+441823.7 


FHO 30 


17.2 


Fader 


Short 0.6 mag fades, typically 





2 days or less, separated by be- 
tween 10 and 60 days. Two 
0.15 mag fades lasting 30 days 
each in mid-2011 and late 2012. 
All fades are superimposed on 
0.4 mag erratic variability. 



Note. — Rmed denotes the median PTF magnitude. Light curves for all these sources are available online from the PTF 
website. 



